How to measure the spreading width for decay of superdeformed nuclei 
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A new expression for the branching ratio for the decay via the El process in the normal-deformed 
band of superdeformed nuclei is given within a simple two-level model. Using this expression, the 
spreading or tunneling width for superdeformed decay can be expressed entirely in terms of 
experimentally known quantities. We show how to determine the tunneling matrix element V from 
the measured value of and a statistical model of the energy levels. The accuracy of the two-level 
approximation is verified by considering the effects of the other normal-deformed states. 

PACS numbers: 21.60.-n, 21.10.Re, 27.70.+q, 27.80.-|-w 



m ; 
o . 
o ■ 

(N ; 

< 



> 

o 

o 
m 
o 



X 



Since the first discovery of superdeformation in ^^^Dy 
Q , one of the principal challenges has been to develop a 
consistent theory regarding the decay-out mechanism of 
the superdeformed (SD) rotational band into the normal- 
deformed (ND) band. Although much experimental 
progress has been made since this first discovery, e.g., 
Lauritsen et al. and references therein, no consis- 
tent theory has been achieved for this decay-out process, 
and, in fact, considerable confusion still exists regarding 
the application of the current theoretical interpretations. 
The purpose of this Letter is to report on new theoretical 
developments, which permit a direct determination of the 
spreading (or tunneling) width for decay out of the SD 
band in nuclei in terms of experimental quantities and, 
thereby, to obtain the magnitude of the tunneling matrix 
element. 

In an earlier publication [i^ , two of us presented a sim- 
ple two-level model to explain the decay out of the SD 
band in nuclei. Employing a retarded Green's function 
approach, we obtained an exact solution for the branch- 
ing ratio, Fm, for decay via the El process in the ND 
band, in terms of the decay widths Fg and F^r in the 
SD and ND potential wells, respectively; V , the tunnel- 
ing matrix element connecting the SD state with the ND 
state, i.e., the state with which it mixes most strongly; 
and A = en — es, the difference between the unperturbed 
energies of these two states (see Fig.^for a graphical rep- 
resentation of these quantities). Our result yielded 
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where f — (F5 + FAr)/2 and V is taken to be positive 
definite without loss of generality. 

We have recently observed that Eq. Q can also be 
rewritten in the form 
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where 
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FIG. 1: Schematic diagram of the two-level problem. V is the 
tunneling matrix element connecting the two states. Electro- 
magnetic decay within each band gives the states their finite 
widths Fiv and Vs. £S and en axe the energies of the two 
levels in the absence of V. 



is the correct expression from Fermi's golden rule for the 
spreading (or tunneling) width 3] . Eq. ^ clearly shows 
that decay into the ND band is a two-step process, and 
that F-^ is a real, physical rate, not a mere theoretical 
construct. Note that 



F^ ^ 27r{V^)/DN, 



(4) 



as was employed by Vigezzi et al. 13 . In Eq. Q, 
is the mean square of the coupling matrix elements con- 
necting the SD and ND states and Dn is the average 
spacing of the ND states. In fact, Eq. gives the av- 
erage spreading width over a fiat distribution in A j^, 
which can deviate drastically from the exact value given 
by Eq. ®. 

Importantly, F-'-, as given by Eq. 0, is a measureable 
quantity, which can be determined from Eq. 



F^ = 



Fn^n^s 



rN-FN{TN+Ts)' 
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F^ is not a free parameter as was stated in Ref. |3- 

We have used Eq. to determine F-^ from the values 
of F/v, F5, and Tn given in Table I of Ref.H and obtain 
values differing radically from those listed in Ref. 0, by 
three to five orders of magnitude (see our Table P). This 
is a direct contradiction of the statement in Ref. that 
"For cases encountered in experiments, their results (i.e.. 
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TABLE I; Tunneling widths extracted from Eq. com- 
pared with the resuhs given in Table I of Ref. 2: (denoted by 

r^'^'). Fm = Pout, rs,rjv, and Dn are the same as those in 
Table I of Ref. Q. The spin values of the decaying states are 
given in parentheses. 



Nucleus 


Fn 


Ts 


Tn 


Dn 


pi 


ri(2) 




— Pout 


(meV) 


(meV) 


(eV) 


(meV) 


(meV) 




0.40 


10.0 


17 


220 


11 


41,000 




0.81 


7.0 


17 


194 


-40 


220,000 


i3*Hg(12) 


0.40 


0.108 


21 


344 


0.072 


560 


i«4Hg(10) 


0.97 


0.046 


20 


493 


1.6 


37,000 



the results of our Ref. |^ agree with those obtained with 
the model of Vigezzi et al." i.e., the model that was 
employed in Ref. 0- 

It should be noted that we obtain a negative value for 
pi i52-|-)y J — 26. This is, of course, physically 



impossible: is a positive definite quantity, 
therefore requires 



Fn < 



N 



■ N 



Eq. © 



(6) 



Experimental results in which the inequality © is vi- 
olated may indicate that the experimentally measured 
value of Fpf is too large or that the SD state is mixing 
with a second ND level, since our result is obtained for 
mixing with only one ND state. At the present time, 
however, uncertainties in the known values of r^r and 
Ts, which are of the order of 100% for Fjv and of the 
order of 10% for Ts, or more, mean that cannot be 
meaningfully determined in cases such as ^^^Dy(I=:26). 

In order to determine the tunneling matrix element V 
from via Eq. we must know A, which generally 
is not experimentally known. We therefore compute the 
expected value for A based on the assumption that the 
states to which the SD state decays in the ND well are 
distributed according to a Gaussian orthogonal ensemble 
(GOE). In the GOE, the probable spacing between levels 
is given by the distribution 



Pis) 



(7) 



where s is the spacing in units of Djv- 

In the absence of tunnehng [V — 0), the energy spectra 
of the ND and SD wells are uncorrelated. Given a spacing 
sDn between the nearest ND levels above and below the 
decaying SD level, then, A is given by the rectangular 
probability distribution 
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Here Q is the Heaviside step function, which simply en- 
sures that A is the nearest neighbor. The probability 
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FIG. 2; Probability distributions for the two ND levels 
bracketing the SD level of interest, from Eqs. @ and 1141 . 
Note that the mean nearest-neighbor spacing is (|Ai|) = 
Dm/ 4:, while the mean spacing of the next-nearest neighbor 
is (IA2I) =37?jv/4. 



distribution for A is then 
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Figure 13 shows V{A)Dn plotted as a function of A/Dn. 
From Eq. (jSJ, it is easy to compute the average detuning 
(|A|) = Dn/^. 

Our ultimate goal is to find the probability distribution 
V{V) for given values of F-^, f and Dn, which in general 
is given by V{V) = 2P(A) |^|. From Eq. ® we can 
obtain |A| as a function of V, 



(10) 



V obviously has a lower bound of Vmin — y ^^^F due to 
the requirement that |A| be real. 

Computing V(V) for the allowed region, we find 



ViV > Vm 



2TV7r 

r^\A\DN 



erfc 



where |A| is given by Eq. (|10|l . The average value of V 
is 



(V) 



Dn 
4 
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Dn 



(12) 



Earlier attempts to consider a statistical theory of SD 
decay-out 0, have been handicapped by their focus on 
average values of and Fn. As was shown in Ref. 0, 
the statistical fluctuations of these quantities are much 
greater than their averages, making such computations 
almost meaningless. This is not the case with the dis- 
tribution of V given by Eq. Ullfl . however, for which the 
fluctuations are of order {V). Furthermore, F-'- and Fn 
are experimentally fixed, whereas the uncertainty in V 
is due to lack of knowledge of the unperturbed energies 
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of the ND and SD states. Eq. 1)11(1 therefore represents 
a central result, since it is the maximal information we 
can have about V without prior knowledge of the shape 
of the potential. 

Thus far we have treated only a two-level model of the 
superdeformed decay-out process. It is reasonable to ask 
what effect the inclusion of more ND states could have. 
As a first step, we include the next-nearest neighbor state 
in the ND well. 

We can again turn to the assumption of the GOE in 
order to find typical values for Ai and A2, the energies 
of the nearest and next-nearest neighbor states, respec- 
tively. 'P(Ai) is simply the P{A) given by Eq. ©. The 
calculation of P(A2) is similar to that for 7^(Ai), but we 
must now concern ourselves with the question of whether 
Ai and A2 have the same or opposite signs. 

In all cases of physical interest, we find that for a given 
spacing s, the correction to the two-state result necessi- 
tated by the inclusion of a third level is larger when the 
two ND levels bracket the SD level. This is quite natural, 
since Eq. ^ requires that the next-nearest ND level in 
this configuration lie on average 40% closer to the decay- 
ing SD state than it would if the nearest two ND levels 
lay on the same side of the SD level. Consideration of the 
case in which Ai and A2 have the same sign, then, only 
decreases the necessary correction. Since our goal is to 
set a reasonable upper bound on this quantity, we assume 
the ND levels lie on opposite sides of the SD level. De- 
cays for which this is not true will in general conform to 
the two-state approximation with even greater accuracy. 

Based on this assumption, we can construct a density 
function for A2 similar to Eq. (jHJ, 



sDn \Dn 2y \ Dm 



(13) 



Together with Eq. Q, this yields a distribution for A2 



^(A2) 



2Dn 



rf(V^|A2 



erf(YI^2|) 



(14) 



This expression for 7'(A2) is illustrated by Fig. [21 Its 
average detuning is (IA2I) = SDn /A. 

Having computed the average values of Ai, A2, and 
V, we are now in a position to begin to see the effect 
of a second ND level. In general, Eq. |^ would suggest 
that the contribution to the total branching ratio of a 
second level is substantially less than that of the nearest 
neighbor. Since Eq. was derived in the context of only 
one ND level, however, we ought to seek a more rigorous 
theory for the three-state branching ratio. In particular, 
we should expect that effects such as competition and 
interference will play a role in the exact result. 

The Hamiltonian for the three-state system can be 
taken to be the sum of two parts, Hq, which represents 
the independant SD and ND wells, and V, which mixes 
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FIG. 3: Plot of ND branching ratios calculated in the two 
i-^N^) ^'^'^ three-state (Pj^^) models with parameters relevant 
to the A ~ 190 mass region. The curves labeled F^l and Pj^^ 
represent the branching ratios for the individual states of the 
three-level model, which sum to the total branching ratio. 
The energy levels were taken to lie at their mean detunings, 
and constant values of Fs /Fat — 10~^ and T/Dn = 10 ~* were 
used. These orders of magnitude are consistent with Table 



the states of the two wells. Hq can be written 

Ho = ^ ^ic\ci + Hem, (15) 

i 

where the sum on i runs over S, Nl, and N2, oi is the ani- 
hilation operator for state i, and Hem contains the cou- 
pling to the electromagnetic field which gives the states 
their nonzero widths. Since they occur by the same de- 
cay process, we assume the widths of the ND states are 
equal, i.e. Fati = rAr2 = Fat. 

Taking y as a perturbation, it is a trivial exercise to 
use Dyson's Equation to construct the retarded Green's 
function of the system. The result, exact to all orders in 
V, is given in the \S), \Nl), \N2) basis by 



G-\E) 



-Vi E-Ai 



-V2 




^V2 







(16) 

where Vi and V2 may be chosen positive without loss 
of generality. In the following, we assume further that 

Vi=V2. 

The branching ratios of the full three-state system can 
now be computed from Parseval's theorem 



dE 
2^ 



|(5|G(i?)|z)p 



(17) 



where i = S,N1,N2. These integrals can be done ana- 
lytically by Cauchy integration, but the results are alge- 
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FIG. 4: Plot of ND branching ratios calculated in the two- and 
three-state models with parameters relevant to the A ~ 150 
mass region. The energy levels were taken to lie at their mean 
detunings, and constant values of Fs — Fjv and T/Dn ~ 10"* 
were used. These orders of magnitude are consistent with 
Table|I] The notation for branching ratios used here is defined 
in the caption of Fig. |3 

braically complicated. It is sufficient for our purposes to 
compute them numerically. 

With Ai and A2 determined by their probability dis- 
tributions, the only remaining parameters in the three- 
state problem are Ts, ^n, V, and D^. In all cases 
of physical interest, however, rsjFjv ^ V,D]y, so we 
can restrict the relevant parameter space by varying only 
the combinations of parameters Ts/^n and V/Dpj sep- 
arately. The corrections to the branching ratios required 
by such a restriction are of order Tn/V. 

Furthermore, the order of F5/F7V is determined by the 
mass region of the nucleus (see Table Pl . In the results 
that follow, variation of this parameter over reasonable 
values does not significantly impact the necessary correc- 
tion to the two-state result. 

Figures 01 and 01 show comparisons of the two- and 
three-state branching ratios for the A « 190 and A « 150 
mass regions, respectively. These figures demonstrate 
that in cases of physical interest, the correction to the 
two-state system due to the presence of a third level is 
relatively small. In the A « 190 region, in particular, 
we find that .9 < F^' /F^' < 1, where the superscripts 
on the branching ratios indicate the number of states in- 
cluded in their calculation. In the A w 150 region, we 
find that 0.7 < F^-'/F^-' < 1. The increased importance 
of additional levels in the ^4 w 150 region arises because 
the typical tunneling matrix elements at decay out are 
significantly larger, which follows from Eq. and the 
relative sizes of F5 and F^v in the two mass regions. Con- 
structive interference between the two ND levels, plainly 



visible in Fig. 01 where Fj^' exhibits a distinct maximum 
for V ~ D]s[/2, can also enhance the importance of ad- 
ditional ND levels. 

In the 190 mass region, however, quantum interference 
effects are seen to be negligable, and this approach can be 
extended to set an approximate upper bound on the total 
error incurred by ignoring all of the ND states except the 
nearest neighbor. It is clear from Eq. that, neglecting 
quantum interference effects, a level's branching ratio is 
approximately proportional to . Assuming that each 
level in the infinite ND band lies at its average detuning, 
we can write the total n-level branching ratio as 

We thus have the result 

00 17(2) 2 
- (2fc _ 1)2 - 8 ■ ^^^> 

In the 190 mass region, then, the expected correction is 
no more than about 23% of Fj^\ 

The three-state results, together with the arguments of 
Eqs. H18() and H19|l . demonstrate that the two-state model 
is sufficient to describe the dominant decay-out process 
of SD nuclei. Within the two-state model, we have shown 
that the decay out of an SD level via El processes in the 
ND band is a two-step process, whose branching ratio 
(0) is expressed in terms of three measureable rates, F5, 
Tn, and TK We have also shown how to determine the 
tunneling matrix element V [Eqs. (|ll|l and (|12|l ] from the 
measured value of F^ and a statistical model of the ND 
band. It is hoped that these results will help clarify the 
nature of the decay-out process in SD nuclei. 
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